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Abstract 

Choice of an appropriate (3+l)-foliation of spacetime or a (2+l)-foliation 
of the Cauchy space, leads often to a substantial simplification of various 
mathematical problems in General Relativity Theory. We propose a new 
method to construct such foliations. For this purpose we define a special 
family of topological two-spheres, which we call "rigid spheres". We prove 
that there is a four-parameter family of rigid spheres in a generic Rieman- 
nian three-manifold (in case of the flat Euclidean three-space these four 
parameters are: 3 coordinates of the center and the radius of the sphere). 
The rigid spheres can be used as building blocks for various ("spherical", 
"bispherical" etc.) foliations of the Cauchy space. This way a supertrans- 
lation ambiguity may be avoided. Generalization to the full 4D case is 
discussed. Our results generalize both the Huang foliations (cf. [3]) and 
the foliations used by us (cf. |8j) in the analysis of the two-body problem. 
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1 Introduction 



The choice of an appropriate coordinate system may sometimes be decisive in 
the analysis of dynamical properties of a physical system. Symmetries of the 
system provide here an important hint. Unfortunately, a generic configuration 
of gravitational field in general relativity theory admits no symmetries. But still 
the nature of the problem suggests very often at least a topological character 
of the coordinate system. As an example consider the quasi-local mass which, 
according to Penrose (cf. [ID]), can be assigned to a topological two-sphere ScE, 
where S is a Riemannian 3D manifold, playing the role of a Cauchy surface for 
Einstein equations (cf. [ID])- Here, topologically spherical coordinates, such 
that S = {r = const.} for a certain radial coordinate r, enable us to simplify 
considerably the formulation of the theory. Recently we have proposed a new 
approach to the two-body problem, where topologically bi-spherical coordinates 
are extensively used to describe the gravitational field surrounding two black holes 
(cf. |8j). In both cases, an appropriate family of two-surfaces, homeomorphic to 
S* 2 , has to be chosen as an important geometric ingredient of the theory. 

A possible construction of such a family consists e.g. in imposing an elliptic 
(cf. |2J) or a parabolic (cf. [5], [6]) partial differential equation for the radial 
coordinate r. An obvious drawback of such an approach consists in the fact that 
we do not control intrinsic properties of the surfaces {r = const.} constructed 
this way. 

In the present paper we propose an alternative approach, based on intrinsic 
geometric properties of special two-spheres. Such objects can be later used as 
building blocks for various foliations. We mimic here the properties of the four- 
parameter family of metric spheres in the Euclidean three-space E 3 , which are 
uniquely characterized by a simple requirement: k = const., where k denotes the 
mean extrinsic curvature. In a generic Riemannian three-manifold, the above 
equation is, however, too restrictive and admits a single solution, related to the 
"center of mass" of the geometry (cf. [4]), whereas in the flat case it admits a 
four-parameter family of solutions. This proves that the condition is not stable 
with respect to small perturbations of the geometry. Moreover, even in the (flat) 
Euclidean space this condition may be replaced by a weaker one, which still selects 
the same family of spheres. We show in the present paper that an analogous, 
weaker than k = const., condition admits a four-parameter family of solutions 
in a generic non-flat case, provided the geometry of the space does not differ 
too much from the Euclidean geometry. Topological two-spheres satisfying this 
condition will be called "rigid spheres". We expect that foliating the Cauchy 
three-surface S with various families of rigid spheres, we will be able to construct 
useful gauge conditions in general relativity theory, similarly as in [8]. 

Our construction is based on the following idea. Take any surface S satisfying 
the rigid sphere condition and consider infinitesimal deformations of S. They may 
be parameterized by sections of the normal bundle T ± S. If we want our condition 
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to admit a four-parameter family of solutions, like in the flat case, its linearization 
must admit a four-parameter family of deformations. This means that we are 
not allowed to constrain the complete information about the mean curvature k: 
four real parameters describing k must be left free. In the flat case these four 
parameters which have to be left free are: the mean value (or the monopole part) 
of k, which is responsible for the size of S, and its dipole part which vanishes due 
to Gauss-Codazzi equations. The dipole part of the deformation is related to the 
group of translations. In fact, possible motions of metric spheres are described by 
the group of Euclidean motions, quotiented by the subgroup of rotations which 
are internal symmetries of every particular sphere S. 

To implement the above idea, an intrinsic, geometric notion of a multipole 
expansion on an arbitrary Riemannian two-surface which is topologically S 2 is 
proposed in Section 2 and a notion of a rigid sphere is then defined. In Sec- 
tion 3 we prove that a four-parameter family of rigid spheres exists in a generic 
Riemannian three-space. 

The present paper is a part of a bigger project, where an eight-parameter 
family of similar "rigid spheres" is proven to exist in a generic four-dimensional 
Lorenzian spacetime. In the present paper we limit ourselves to the case of a 
Riemannian three-manifold, but our construction can be generalized to the entire 
pseudo- Riemannian spacetime M, instead of the Riemannian Cauchy three-space 
E C M. The idea of this extension is to mimic the case of the flat Minkowski 
space, where all possible round spheres, embedded in all possible flat subspaces 
E of M, form an eight-parameter family. All of them can be obtained from a 
single one by the action of the product of the one-parameter group of dilations 
(changing the size of S) and the ten-parameter Poincare group, quotiented by the 
three-parameter rotation group. In this case, it is necessary to take into account 
not only the external curvature of S, but also its torsion. The rigid spheres 
obtained this way form an eight-parameter family and may be used to construct 
useful coordinate systems not only on a given Cauchy surface S, but also in 
the entire spacetime. The main advantage of such a construction consists in 
its rigidity at infinity: it eliminates supertranslations and reduces the symmetry 
group of the Scri, otherwise infinite dimensional, to the finite-dimensional one. 
The corresponding result will be presented in a subsequent paper. 

2 Equilibrated spherical coordinates 
2.1 Conformally spherical coordinates 

Let S be a differential two-manifold, diffeomorphic to the two-sphere S 2 C M 3 
and equipped with a (sufficiently smooth) metric g. Coordinates (t?, ip) = (x A ), 
A = 1,2, defined on a dense subset of S\£, where £ is topologically a line interval, 
will be called conformally spherical coordinates if they have the same range of 
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values as the standard spherical coordinates on S 2 C M. 3 and, moreover, if the 
corresponding metric tensor g^B is conformally equivalent to the standard round 
metric on S 2 , i.e. the following formula holds: 



9ab = V> ■ °ab , 
where ip is a (sufficiently smooth) function on S and 



1 

sin 2 tf 



(1) 



(2) 



Remark 1. Conformally spherical coordinates always exist (cf. [9]). It is easy to 
see that there is always a six-parameter freedom in the choice of such coordinates. 
More precisely: if (■#, (p) are conformally spherical coordinates then (i?, (p) are 
also conformally spherical if and only if they may be obtained from tp) via a 
conformal transformation of S 2 C I 3 . 

Example 1. A "proper" conformal transformation, i. e. which is not a rotation: 
Let n G S and r > be a positive number. Using appropriate rotation, choose 
conformally spherical coordinates ip) in such a way that n is a north pole, 
i.e. the coordinate vanishes at n. Define 



where 



or, equivalently, 



F n J0,<p) = (■&,$) 



{} := 2 arctan r • tan — 

2 





tan — = r ■ tan ■ 



ip:=<p 



2 2 

For the fixed point n these transformations form a one-parameter group: 



F o F = F 



generated by the vector field: 



d. 
df 



t=l 



dt 



t=i 



2 arctan I t ■ tan 



d_ 



sin0 



d_ 



(3) 
(4) 

(5) 
(6) 

(7) 



In particular, F Uj i = I (the identity map) for every n. Moreover, equation (J3]) 
implies the following identity: 



F-nr — F n 1 . 

' T 

Using 01]) and (j3J) we may easily derive the following formula: 



d0 ~ 
d0 = d0 
d0 



1 + tan 2 | ~ 

r dd . 

r 2 + tan 2 | 



(8) 



(9) 
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Similarly, we may prove: 

sintf ~ 2tan| 
sin v = ~ sin v = 



sini? 



1 + tan 2 f 



2±tanf 

T Z 



2 i? 

1 + tan 2 t 



1 + 4 tan 2 | 2 tan | 



1 + tan 2 | ~ 
2 tan | 

1 + tan 2 f 
sin = r 4 sin $ 

r 2 + tan 2 f 



As a conclusion we obtain: 

(ch?) 2 + sin 2 (oV) 2 = h 2 [(dtf) 2 + sin 2 (dy?) 2 ] , 



where 



h 



1 + tan 2 f 
r 2 + tan 2 f 



(10) 
(11) 

(12) 
(13) 



which proves the conformal character of the transformation. Indeed, we have 
g AB dx A dx B = $ [(dtf) 2 + sin 2 tf(dy?) 2 ] = ^h 2 \{d$) 2 + sin 2 i?(d^) 2 l . (14) 

Hence, </?) are conformally spherical coordinates if <p) were. 

2.2 A barycenter of a conformally spherical system 

Given a system of conformally spherical coordinates on S, consider the corre- 
sponding three functions: 



x 

y 

z 

We have, therefore, a mapping D 

D(0,v>) = 

The following vector 

X 



= sin $ cosy? , 
= sin^sin^ , 
= COS 1 !? . 

]0,7r[x]0,27r[H> 

D\$,<p) 
D 2 (^<p) 



(15) 
(16) 
(17) 




given by: 




(19) 



where by < / > we denote the average (mean value) of the function / on S, i.e. 
the number 

j. J s fy/det g d 2 x 
f s \/det g d 2 x ' 

will be called a "barycenter" of the system if) on S. Of course, we have 
11X11 < 1. 
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Example 2. Consider the proper conformal transformation 
new barycenter 

< x > 



X 



where 



x 

y 

z 



< l > 

< z > 



sin i? cos f , 
sini9sin<^ , 
cos?? . 



and calculate the 



(21) 



The trigonometric identity: 



cos $ 



1 - tan 2 f 
1 + tan 2 | 



implies: 



tan 







COS0 



1 + COS 1 + -2 



Hence, the formula (j^J) implies: 



or, equivalently, 



2 



(22) 
(23) 

(24) 
(25) 

(26) 

(27) 

1 + Z + T 2(!_ z y ■ ( 28 ) 

To calculate mean value of the functions (12 7p . (128)) and (|25p we do not need to 
pass to new coordinates (0, ip), but we may use, as well, old coordinates (0, <p). 
But we see that for r — >■ we have x — > 0, y — >■ 0, z — > 1. The Lebesgue theorem 
implies, therefore, that for r — >• we have 



1 + z 1 + z 

1 + Z - T 2 {\ - Z) 

" 1 + z + r 2 (l - 2) 
Moreover, the formula ffTTjl and its inverse: 

1 + tan 2 f 1 + £f 
sin0 = r — -r sint? = r -? — sin0 



2r sin 



r 2 tan 2 ^ 



r 



2 H l + 2 + r 2 (l-2) ' 

J- 1 -2 



imply: 



2r 



l + 2 + T 2 (l 

2r 



a; 



V ■ 





( < £ > ^ 






H 


<y> 


HI 


< > 




\ <z> J 







n . 



(29) 
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2.3 Equilibrated spherical systems 

Definition 1. Conformally spherical coordinate system ($,</?) is called equili- 
brated, if its barycenter vanishes: X = £ I 3 . 

Remark 2. If there are two equilibrated spherical systems on S then they are 
related by a rotation. 

Theorem 1. Each metric tensor on S admits a unique (up to rotations) equili- 
brated spherical system. 

Proof. Given a metric tensor g on S, choose first any system of conformally 
spherical coordinates tp) on S and consider the corresponding identification of 
its points with the points of S 2 = dK(0, 1) C M 3 . Consider now the mapping 

M 3 D #(0,1) 3 N -»■ F(N) G X (0,1) C M 3 , (30) 

given by the following formula 

J-(N) := X n , T , (31) 

where the latter is the barycenter of the coordinates (p) obtained from </j) 
by the proper conformal transformation Q with 

n := PI (32) 

and 

r := 1 - ||N|| . (33) 

Obviously, F is continuous. Moreover, due to it reduces to the identical 

mapping when restricted to the boundary S 2 = dK(0, 1) C K(0, 1) (i.e. we have 
J-'(N) = N for ||N|| = 1). This implies that there must be a point No which 
solves equation J 7 (No) = 0. This completes the existence proof. To prove the 
uniqueness, let us suppose that there is another solution: .F(Ni) = 0. Consider 
now the conformal transformation F nijT1 o F^ TQ . Since the proper conformal 
transformations do not form any subgroup of the group of all conformal transfor- 
mations, we cannot assume that it is again a proper transformation. But it may 
be decomposed into a product of rotations and a proper conformal transforma- 
tion: 

Fn 1 ,r 1 oF n ^ T0 = O 1 oF m , T oO , (34) 

where 0\ and Oo are rotations. Denote by ($ , ip ) the spherical coordinates ob- 
tained from ($, if) by the transformation F no<TO and then rotation Oq 1 . Similarly, 
denote by ^i) the ones obtained from ($,</?) by F Ul)Tl and then by rotation 
Oi . Because a rotation does affect equilibration of coordinates, both systems 
(i^Oj^o) an d are equilibrated. But the latter may be obtained from the 
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former by a proper conformal transformation -F mT . We shall prove that this is 
impossible unless r = 1 or, equivalently, transformation F mr is trivial (identical). 
For this purpose consider, for each value of r, the linear function z T . Without 
any loss of generality we may assume that 



m= (35) 




(if this is not the case, it is sufficient to perform an appropriate rotation of 
coordinates). Formula (125 p implies the following relation: 



1 + Z -T 2 (l- Zq) 
ZT 1 + Z + T 2 (1-Z ) 

Hence 



(36) 



d_ _ -At {I- z 2 
dt Zr ~ [i + Zo + r 2 (1 _ Zq)] 

and it vanishes only at a single point z = 1. Consequently, its mean value: 



< 1* > (38) 

is strictly negative. This implies that starting from r = 1 (which corresponds to 
the identity mapping -F m ,i) and moving towards the actual value r < 1, the "z"- 
component of the vector X m r is strictly increasing. It vanishes at the beginning 
because ($0)¥>o) is equilibrated. Hence, it must be strictly positive at the end. 
This means that the final system ($i,<^i) cannot be equilibrated unless r = 1 
and, therefore, both systems coincide. □ 

Different equilibrated spherical systems of coordinates form, therefore, a three- 
dimensional family. They can be parameterized by the position of a fixed point 
n G S (north pole) and the geographic longitude of a fixed point m G S (Green- 
wich). More precisely: given two points n, m G S, n ^ m, there is a unique 
equilibrated spherical system (p) of coordinates on S, such that # vanishes at 
n and if vanishes at m. 

Combining these observations with classical results (cf. |9]) we obtain the 
following 

Theorem 2. Let S be a differential two-manifold, diffeomorphic to the two-sphere 
S 2 C M 3 and equipped with a metric g of class C^ k ' a \ For every pair n, m G S, 
n/m, there is a unique equilibrated spherical system ip) of coordinates on S , 
such that $ vanishes at n and tp vanishes at m and the metric components qab 
are of the same class C^ k ' a ^ . 
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Here, is a Holder space C k ' a (S 2 ), defined for 1 ^ k 6 N and < a < 1. 

The space consists of those functions on S 2 which have continuous derivatives up 
to order k and such that the A;-th partial derivatives are Holder continuous with 
exponent a. This is a locally convex topological vector space. 
The Holder coefficient of a function / is defined as folows: 

\J\C0," = SUp _ . 

x,y&S 2 ,x^y \ x U\ 

The function / is said to be (uniformly) Holder continuous with exponent a if 
|/|co,a is finite. In this case the Holder coefficient can be used as a seminorm. 

The Holder space C k,a (S 2 ) is composed of functions whose derivatives up to 
order k are bounded and the derivatives of the order k are Holder continuous. It 
is a Banach space equipped with the norm 

||/||c^ = ||/||^+max|^/| c o,«, 

\P\=k 

where (3 ranges over multi-indices and 

II /lie* = max sup \D p f(x) \ . 

\P\<k x€ s 2 



2.4 Rigid spheres in a Riemannian three-space 

Given a manifold S equipped with a metric tensor g, there is a three-dimensional 
space of "linear functions" uniquely defined on S as linear combinations of func- 
tions (|T5T - [T7|) . calculated in any equilibrated spherical system of coordinates 
("&,(p). We denote this space by Ai 3 . By Ai 4 we denote the space spanned 
by .M 3 and the constant functions on S. Linear functions (IToTfTTj) on S are 
eigenfunctions of the Laplace operatoi0 A CT , with the eigenvalue equal to —2, i.e. 
ApX 1 = —2X\ where we denote x = X 1 , y = X 2 , z = X 3 . Let us denote by 
da := sin $ d?9 d(p the measure associated with the metric (Jab- 

Definition 2. Let / G L 2 (S,da). The projection of / onto the subspace of 
constant functions: 

P m {f) := ^ ! f da (39) 
J s 

will be called the monopole part of /, whereas the projection onto 
M 3 = span{X 1 ,X 2 ,X 3 }: 

1 By A CT we denote the usual Laplace operator for the unit-sphere metric ©. 
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will be called the dipole part of /. In addition, we set 

M 4 := span{l} © M 3 = span{l, X 1 , X 2 , X 3 , } (41) 

and Pmd(f) '■= Pm(f) + Pd(f) £ -M 4 denotes the mono-dipole part of /. 

The above structure enables us to define the multipole decomposition of the 
functions defined on a topological sphere S in terms of eigenspaces of the Laplace 
operator associated with the metric (Tab- If h is a function on S then by h m : = 
P m (h) we denote its monopole (constant) part, by h d := Pd(h) — the dipole part 
(projection to the eigenspace with 1 = 1). By h w := (I — P m d)(h) = h — h m — h d 
we denote the "wave", or mono-dipole-free, part of h, h dw := (I — P rn )(h) = 
h-h m = h d + /i w , and finally h md := P md (h) = h m + h d . 

Remark 3. Mutually orthogonal projectors P md and P w := (I — P md ) are, of 
course, continuous, when considered as operators in the Hilbert space L 2 (S, da). 
For our purposes we have to consider them as operators in the Banach space 
C^- k ' a \ Here, no "orthogonality" is defined. Nevertheless, both operators are 
again continuous projectors. They define an isomorphism: 

s~i(k,a) ~ /nf(fc,a) Wfc.a) 
° — °md X U w i 

where C^f = P md {C^) = M A and C ( w k ' a) = P w {C (k ' a) ). Hence, a function / G 
(j{k,a) j g un jq Ue iy characterized by its mono-dipole part / md and the remaining 
"wave" part / w , i.e. we have: / = (/ md ,/ w ). 

Definition 3. Let S be a Riemannian three-manifold and let S C £ be a sub- 
manifold homeomorphic with S 2 C M 3 . We say that S is a rigid sphere if its 
mean extrinsic curvature k satisfies k G Ai 4 , i.e. if the following equation holds: 

k w = . (42) 



2.5 The 4-D spacetime case - an outline 

Definition of a rigid sphere in a Lorenzian four-manifold is more complicated: to 
control "rigidity" of a sphere, we must take into account more geometry. For this 
purpose we consider the extrinsic curvature vector of S: k a = k AB g AB , where k AB 
denotes the external curvature tensor of S (here, a, b are indices corresponding to 
the subspace orthogonal to S whereas A, B label coordinates on S) . Moreover, 
we consider its torsion: 

t A = (m|V A n) , (43) 

where 

n := H ' (44) 

||k|| = a/ k a g ab k b , and m is a vector orthogonal to both k and S. 



10 



Definition 4. Let M be a Lorenzian four-manifold (a generic curved spacetime) 
and let S C M be a spacelike submanifold homeomorphic with S 2 C K 3 . We say 
that S is a rigid sphere if k = {k a ) is spacelike and the following two conditions 
are satisfied: 

||k|| G M 4 , (45) 
V A t A e A^ 3 . (46) 

In this paper we limit ourselves to the purely Riemannian 3D-setting. The 
general, pseudo-riemannian case will be analyzed in a subsequent paper. 



Example 3. Rigid spheres in a four-dimensional Minkowski spacetime and in 
Euclidean three-space. 

Let Mq be the flat Minkowski spacetime, i.e. the space M 4 parameterized by 
the Lorentzian coordinates ( ) and equipped with the metric 

V — (Vap) — diag(— 1, 1, 1, 1) (greek indices run always from to 3). 
Consider in M a round sphere, i.e. the two-dimensional submanifold defined by 



Sr,R:={xeR A 



x° 



T , J>*) 2 = i? 2 } 

i=i 



where the time T e K. and the sphere's radius R > are fixed. It may be easily 
verified that the submanifold fulfills the following conditions: 



Vk a g ab k b = I G M 4 , (47) 
Va^ 4 = G 7W 3 , (48) 

hence each round sphere St,r in Minkowski spacetime M is a rigid sphere. Using 
Poincare symmetry group of M , it is easy to check that there is an 8-parameter 
family of such spheres. Indeed, fixing the value of R, a 7-parameter family remains 
left. All of them may be obtained from a single sphere, say So t R, by the action 
of the 10-parameter Poincare group. Because the three-parameter subgroup of 
rotations corresponds to internal symmetries of So,r, we are left with 7 parameters 
only. The parameter R corresponds to the dilation group. Hence, we have 10 — 

3 + 1 = 8. 

In Euclidean three-space (represented by a slice {x° = 0} in M ) the family 
of rigid spheres reduces to four-parameter family of such spheres, where 4 = 3 + 1 
- three translations plus dilation (or similarity transformations minus rotations 

4 = 7 — 3). Each round sphere in Euclidean three-space is a rigid sphere because 
its mean extrinsic curvature k = — -| G M. 4 . 
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3 Existence of rigid spheres in a Riemannian space 



Let X be a three-dimensional Riemannian manifold. Let S C X be a two-manifold 
diffeomorphic to the unit sphere S 2 C M 3 . We consider the following problems: 
1) Can we deform S in such a way that the resulting submanifold becomes a rigid 
sphere? 2) How many of such deformations exist in a vicinity of 57 

To parameterize these deformations we introduce in a neighbourhood of S 
a Gaussian system of coordinates (u,x A ). Here, by (x A ), A = 1,2, we denote 
any coordinate system on S, whereas u is the arc-length parameter along the 
u {x A = const.}" geodesies starting orthogonally from S. The three- metric takes, 
therefore, the form 

g = du 2 + g AB (u, x A ) dx A dx B . (49) 

Suppose, moreover, that coordinates (x A ) = (■#, (f) are conformal and equili- 
brated on S. This means that we have 

c/ab dx A dx B = i) ■ {a AB dx A dx B ) = tP ■ ( dtf 2 + sin 2 d^ 2 ) , (50) 

where 

°9ab ■= 9ab(0,x a ) (51) 

is the induced two-metric on S, a is the "round" two-metric on the Euclidean 
unit sphere: 

a AB dx A dx B = dtf 2 + sin 2 dip 2 , (52) 

and the function ip is dipole-free (ip d = 0). Second fundamental form of S is 
given by: 

k A B = ~-9ab,u ■ (53) 

Its trace does not need to belong to the space M. 4 of mono-dipole-like functions, 
i.e. the surface S does not need to be a rigid sphere. We are looking for such 
deformations of S, for which the resulting surface fulfills already the rigidity 
condition. 

Any deformation of S which is sufficiently small may be uniquely parame- 
terized by a function r = r(x A ), such that the deformed surface S T is given 
by: 

S t = {(u,x a )\u = t(x a )}. (54) 
The surface S T carries the induced metric: 

g\s T = [dr(x A )] 2 + g AB (r(x c ),x c ) dx A dx B =: g AB (x c ) dx A dx B , (55) 

where 

g AB (x c ) = (d A r)(d B r) + g AB (r(x c ),x c ) . (56) 

Here, we use the same coordinate system (x A ), which was previously used for S. 
However, these coordinates do not need to be neither conformally spherical nor 
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equilibrated. To verify that the deformation r was successful, i.e. that S T is a 
rigid sphere, we have to pass to an equilibrated system of spherical coordinates, 
say x A , on S T . To make this choice unique, we use the north pole: n := {# = 0}, 
and the "Gulf of Guinea": m :={??= ~; <p = 0} to get rid of the rotation non- 
uniqueness (cf. Theorem [2]). This way we obtain an equilibrated version 7}ab of 
the metric (1561) . Finally, we calculate the extrinsic curvature k and check whether 
or not its wave part k w (S T ) satisfies condition k w (S T ) = 0. 

The idea of our paper may, therefore, be sketched as follows. We begin with 
a metric (H9|) which is of the class C^ k,a \ The above construction defines a 
continuous mapping: 

C^ 1 '^ x Clf +1 ' a) 3 (r md , t w ) = t — y F{r) := k™ G C%~ 1 ^ . (57) 

Indeed, the resulting metric in a neighbourhood of S T is obtained from g and 
the first derivatives of r. The function r being of the class C^ k+l,a \ the metric 
obtained this way is again of the class C^ k,a \ Due to Theorem [2], its equilibrated 
version Tjab is again of the same class. Finally, the extrinsic curvature k is 
obtained using first derivatives of this metric. Hence, the result is of the class 
(j(k-i,a) an j £ ne entire procedure is continuous. 

Now, rigid spheres are those, which satisfy equation: 

F(t) = . (58) 

We are going to prove that, for a generic metric g, the above equation defines an 
implicit function: 

M 4 = Cf+ 1,a) 9 r md — ► H(r md ) e Clf +1 ' Q) , (59) 

such that 

F(r md ,#(r md )) =0 , (60) 

or, equivalently, that S^md^^mdi) is a rigid sphere. The main result of our paper 
follows as a corollary: 

Theorem 3. Generically, there exists a four-parameter family of rigid spheres in 
a neighbourhood of a given two-sphere S C S, corresponding to the four-parameter 
family of mono-dipole functions r md on S. 



3.1 Infinitesimal deformations of spheres 

To prove existence of the implicit function (160]) it is sufficient to show that, given 
a mono-dipole deformation r md , the partial derivative of F with respect to the 
"wave-like" deformation r w is an isomorphism of Cw 1 ' ^ onto Cw~ l, °^ ■ 

For this purpose, we analyze the infinitesimal, linear version of the construc- 
tion discussed above. Consider, therefore, a transversal deformation r = t(x a ) 
of S C £ and a small deformation parameter e: 

S ~» S t = {{u,x a )\u = et{x a )} . (61) 
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Under such transformation the induced metric changes in the following way: 

9ab - 9ab = -2eT°k AB + 0(e 2 ) . (62) 

Even if the initial system of coordinates was equilibrated, the transformed metric 
does not need to be conformally spherical. The non-sphericality of the metric 
must be, therefore, compensated by a change of coordinates. Its infinitesimal 
version is described by a tangential (with respect to S) deformation 

x A = x A - ei A . (63) 

Under such coordinate transformation the metric changes as follows: 

9ab = 9ab ~ £ 2£ g 9ab , (64) 

where the last term represents the Lie derivative of the metric g AB with respect 
to the vector field "— on S. But, according to fl62|) . the difference between 
gAB and qab is already of the first order in s. Hence, if we replace it by the Lie 
derivative of the metric gAB, the error will be of the second order in s. Using the 
Killing formula for the Lie derivative of the metric we finally obtain: 

gAB = gAB + 2ef(A||B) + 0(e 2 ) , (65) 

and the covariant derivative iu is taken with respect to the original metric g A s- 
Hence, we have: 

gAB - gAB = -2erk AB + 2e^ a \\b) + 0(e 2 ) . (66) 

Let us decompose the above equation into the trace and the trace-free parts, 
calculated with respect to g A B (we omit the terms of order e 2 and higher): 

gAB - gAB = (z£, C \\c - £T°kj g AB - 2erk AB + 2e (^(a\\b) ~ '^ C \\c9AB S j , (67) 
where 

kab ■= k AB - -kg AB (68) 

is the traceless part of k AB . We want ~g A B to be conformally spherical, i.e. *g AB = 
a ■ g AB . This implies: 

(l-erk- a + e^ C \\cj g AB - 2erk AB + 2e£( A \\ B ) - ^ C \\c9ab = . (69) 
The trace part of this equation defines uniquely the value of a: 

a = l-erk + ef llc , (70) 
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whereas the trace-free part reduces to: 

€a\\b + €b\\a — £ \\c°9ab — 2rk AB . (71) 

It is convenient to rewrite equation (17TI) in terms of the "round" unit-sphere 
geometry <Jab- For this purpose we use the following conventions: components 
of a vector (i.e. an object having upper indices) are the same in both geometries 
oab and cjab = ^^ab- Components of a co- vector (lowered indices) are denoted 
as follows: 

& = <tab£ B , U=°9abZ b = tI>vab = tI>F a . (72) 
The covariant derivative with respect to gab will be denoted by u a, e.g. £^, B . 
Equation (ITT]) can be easily rewritten as: 

CauB + CbuA - £ llCO~AB = ^° K AB ■ (73) 

The left-hand side of this equation defines a mapping from the space of vector 
fields on the unit sphere to the space of trace-free rank 2 tensor fields. The 
kernel of this mapping consists of the dipole The "Fredholm alternative" 

argument shows that the operator on the left-hand side defines an isomorphism 
between the space of dipole-free vector fields on the unit sphere and the space of 
trace-free rank 2 tensor fields (see also |7]). This isomorphism (in metric a) will 
be denoted by i 12 . Hence, the wave part of £ A is implied uniquely by equation 
(I73"|) (see Appendix): 

C A = i£ fy'KAa} , (74) 

whereas the dipole part of £ A , i.e. the field £ dA , remains arbitrary. 

The above choice of the wave-like component of the tangential deformation 
( wyt guarantees that the new coordinate system x A is conformally spherical. We 
would like it to be also: 1) equlibrated and 2) satisfying conditions related to 
the two fixed points n and m. These conditions mean that the field £ has to 
vanish at the north pole n and that its 99-component vanishes at m. The above 
3 + 3 = 6 conditions fix uniquely the total dipole-part of the tangential (to S) 
deformation £ . This way the continuous mapping which assigns uniquely the 
tangential deformation £ A to its transversal component r has been defined. 

3.2 The infinitesimal change of the extrinsic curvature 

Now, we are going to calculate the infinitesimal change of the wave part /c w of the 
mean curvatur^l, i.e. derivative of the mapping (l57j) with respect to the "wave- 

2 A vector field on the sphere may be uniquely decomposed into the sum of a gradient and a 
co-gradient. These two components are represented by the corresponding two scalar functions: 
the divergence and the curl. The multipole expansion of a vector field is uniquely defined by 
the multipole expansion of these two functions. 

3 First variations of the total mean curvature k is known in the literature as the second 
variations of area, cf. e.g. p]. See also discussion in the Appendix. 
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like" deformation r w . We have k = g AB kAB, where g AB denotes the inverse 
of the two-metric qab (whereas g AB denotes the corresponding components of 
the inverse three-metric.) The simplest way to calculate this change is to use a 
coordinate system (u,x A ), adapted to the deformed surface: 

to — u — er(x A ) , i.e. S T = {u = 0} , (75) 

and the formula: 

kAB = ^ W AB • (76) 

The three-metric g takes now the following form: 

g= du 2 + 2st a du dx A + g A s dx A dx B + 0(e 2 ) . (77) 
This implies g^ = 1 + 0(e 2 ) and, consequently, 

k A B = T^ab + g^TcAB + 0(e 2 ) = - (g wA \\ B + g uBU - Qab,w) + 0{e 2 ) , (78) 

where we treat the "shift vector" g^A — ^ T ,A as a covector field on S T . The first two 
terms combine to £T\\ab, whereas the last one: gAB,ui{S T ) can be approximated 
by the quantity gAB,u{S) — ~2kAB plus the derivative of this object. Finally, we 
have 

kAB = k A B + £rk A B,u + st\\ A b + 0(e 2 ) . (79) 

Since the derivative g A B,u> of the metric gAB is described by — 2k ab, the derivative 
of its inverse g is described by +2k . Hence, we have: 

~AB _ <>AB = ^AB + 0(£ 2) ^ (g0) 

and, consequently: 

k = g AB k AB = k + erd u k + et u a + 0(e 2 ) . (81) 

The quantity rd u k + t^ a a describes already the second variation of area (see 
Appendix), i.e. the derivative V T k. However, to calculate the derivative of the 
mapping (j57j) . we have to select its wave part k w . For this purpose we have to pass 
to the conformally spherical, equilibrated coordinates x A , given by formula (163]) . 
Infinitesimal change of the scalar function k with respect to this deformation is 
given by formula: 

k — k — ei A k A + 0(e 2 ) . 

Hence, we get: 

\(k-k}= rd u k + t^ A a - i A k A + 0(e) , (82) 
or, equivalenty (cf. Appendix), 

I (k - fc) = t(R\ + k AB k AB ) + tW A a - i A k, A + 0(e) . (83) 
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3.3 Proof of the Theorem g] 

The last formula gives, finally, the value of the derivative of the mapping (157|) . 
When restricted to the subspace of wave (i.e. mono-dipole-free) deformations, it 
gives us: 

C (k+l,a) 3T ^^ r (^ R u u + k AB kAB j + t \\A a _ £ C (fc-l,a) (g4) 

The above linear operator is, obviously, continuous. In particular, the vector field 
£ A is given by formula (ITlj) . together with the accompanying vanishing conditions 
at n and m. 

If the space E is flat (Euclidean) and S is a standard (rigid) sphere of radius 
r, then we have: 

9ab = r 2 a AB ; k AB = -ra AB ; k A = ; R u u = . (85) 
Hence, the above operator reduces to: 

r w ^ [r(R u u + k AB k AB )+T\\ A A -i A k A \ = y 2 [(A CT + 2)r] w 

= ^(A CT + 2)(r w ), (86) 

which is obviously an invertible mapping from C { w k+1 ' a) to C ( *- l ' a) . But the map- 
ping ( I86p depends in a continuous way upon the geometry (metric and curvature) 
of S. This implies that it remains invertible for sufficiently small deformations of 
the geometry. 

Suppose now, that S is asymptotically flat, i.e. there is a coordinate chart 
(x h ) which covers the exterior of a compact domain D C S, such that 

9kl = $kl + hki , 

where h vanishes sufficiently fast at infinity. This means that the above statement 
remains true for "coordinate spheres" defined as follows: 

3 



S SojR := {x e R 3 | J2( xt - x o) 2 = r2 } 

i=l 



if the radius R is sufficiently big. We conclude that £\D admits a four-parameter 
family of rigid spheres, similarly as in the case of the flat metric. 



4 Conclusions 

The main technical ingredient of this paper is the intrinsic, coordinate invariant 
definition of the "multipole expansion" of a function defined on a Riemannian two- 
manifold, diffeomorphic with S 2 . This enables us to select a finite-dimensional 
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family of "rigid spheres". The dipole part k d of the curvature parameterizes 
the position of the center of such a sphere with respect to the center of mass. 
In particular, A; d = corresponds to the spheres, which are centered at the 
center of mass. Properties of such a foliation have been analyzed in [4j. General 
topologically spherical coordinates, having property that surfaces {r = const.} 
are rigid, do not admit supertranslations ambiguity at space infinity. This way 
symmetries of the "tangent space at infinity" reduce to a finite-dimensional one. 
The 4D version of our results, valid for a generic four-dimensional Lorenzian 
spacetime, which will be presented in the subsequent paper, will do the same job 
for the symmetry group of the Scri. 
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A Appendix 

A.l The dipole part of traceless symmetric part 

The kernel of the mapping 

defined by the left-hand side of the formula (I73p consists of the dipole fields. This 
is a simple consequence of the following observations. 

• In case of the unit sphere the Hodge decomposition £ = da + 5(3 + h of the 
covector £ on a compact manifold does not contain the harmonic part, i.e. 
harmonic one- form h vanishes ( dh = = 5h implies h = 0). The topology 
of the unit sphere (triviality of the corresponding cohomology class) cancels 
the harmonic part and we can always represent £ as follows 

£4 = a A + e A B (3, B , (87) 

where functions a and (3 are defined up to a constant but their gradients 
are unique. 

• The purely dipole covector £ simply means that the potentials a and (3 are 
purely dipole functions: a = a,iX\ (3 = biX\ where a^, bi are real constants. 

• Direct computation for dipole functions X 1 enables one to check the follow- 
ing identity: X^ AB = —X % oab, hence for any dipole function a 

a U AB = -aa A B ■ (88) 
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• Formulae flBZD and (|SSJ) give 

hence the traceless symmetric part of £aub vanishes. 

A. 2 The isomorphism between covector fields and symmet- 
ric traceless tensors on (S* 2 , gab) 

Let us consider the following diagram: 



V k+2 © Vk+2 ) 



yl JHi 
v k+l } 



V k+2 



Fl 



T/0 

V k+2 J 



Ha, 



v k 2 



v k 2 



121, 



HO, 



V k-2 



VP 

V k-2 



V k-2 



Fl 



V k-2 



where the mappings and the spaces are defined as follows: 

ioi(f,g) = fnA + £a B 9ub , 



c 



vaub + v Bu a - &ABV lie 



Ho 



*2i(x) = Xa 
(v) = (i 



UB 



vaub) , 



F Kf> 9) = (9, f), v A = e A B v B , Xab = e A C XcB 



k,a 



k.a 



V k ° - scalars on S 2 belonging to Holder space C k,a 
V k l - covectors on S 2 belonging to Holder space C 
V k 2 - symmetric traceless tensors on S 2 belonging to Holder space C 
Denote by A a the Laplace operator on S 2 and by SH the space of spherical 
harmonics of degree I, (/ G SH 1 -<=>- A a f = —1(1 + l)f). The following equality 

ho ° «2i 112 ° ioi = A CT (A CT + 2) 

shows that if we restrict ourselves to the spaces V := V° [SH° © SH 1 ] = 

{I-P md )V° (A a (A a + 2)V° = V°) andF 1 = V 1 Q [iq^SH 1 )] ((A ct + I)V l = V 1 ) 
then all the mappings in the above diagram become isomorphisms. 



A. 2.1 Integral operators, generalized Green's functions 

Solution of the Helmholtz equation on a unit sphere S 2 : 

[A CT + 1(1 + 1)] *,(n) = $(n) , n e S 2 
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is given (see e.g. [TT]) in terms of the generalized Green's function Gi as follows: 



Gt(n,ri)$(ri)d 2 ri 



(89) 



S 2 



Here n = D($, ip) given by ffT8|) and d 2 n = da. The solution \I/;(n) is automati- 
cally orthogonal to the space SH l (the kernel of Helmholtz operator A a + 1(1 + 1)) 
because Green's function is orthogonal to this space. In our case we need to write 
the inverse of the operator A (T (A cr + 2) as a double integral with the correspond- 
ing kernels Gi for / = and 1 = 1. More precisely, the solution g of the equation 
A (T (A - + 2)g = f (with P m df = 0) is given in the following form: 



g(n) 



P w / G (r*,n" 
Is 2 



Us 2 

Go(n, n 
s 2 



Gi(n",n')/(n')dV / G^m, n')f(n') dV d 2 m 

47r Js 2 xS 2 



Gi(n",n')/(n')dV 



Us 2 



dV , 



d z n 



(90) 



where the projection operator P w provides orthogonalitj0 of g to the space SH°( 
SH l . The generalized Green's function written in a standard form: 



Gl{n,n>)= E J m ^ )Yim , (n '\^ y.n.eSH' 



i=0,i=£l m=—i 

can be simplified as follows (cf. |11|): 
1 



l(l + l)-i(i + l) 



Gi(n,n) = —Pi(n- n') 

47T 



1 — n • n' 
In h Q 



9tt ^ 



2i + 1 



2/ + 



27r^(/-*)(Z + 2 + l) 



2z + 1 



Pi(n • n) , 



i=0 



(Z-z)(J + z + l) ' 



Yi m - spherical harmonics (orthonormal basis in SH l ), n ■ n' G [—1, 1] is a scalar 



product of unit vectors in R 3 and Pi(x) 
mial. 



2'Z! 



— l)^ is the Legendre polyno- 



A.3 Second variation of area 

The Gaussian coordinates (H9|) and the definition of the Riemann tensor gives 

R U AuB = k A B,u + k A C k B c ■ (91) 



4 The above integral operators do not mix the wave part with the mono-dipole part of a 
function. This means that P w f = f implies P W (G\ * f) = G\ * f and P w (Gq * f) = Gq * f ■ 
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This leads to 



k AB = k AB + et(R u AuB - k A c k BC ) + £Tu AB + 0(e 



(92) 



Taking the trace (and using 



we obtain: 



k = k + et(R u u + k AB k AB ) + er" A A + 0(e 



(93) 

The formulae (|8Tj) and fl93|) are equivalent because of the Gauss-Codazzi equa- 
tions: 

2d u k = R(g kl ) - R(g AB ) + k AB k AB + k 2 , (94) 

2R\ = R(g kl ) - R(g AB ) - k AB k AB + k 2 , (95) 

where R(gu) and R(g AB ) are scalar of curvatures of the three-metric g k i and the 
two- metric g AB , respectively. Obviously flSTT) . f[95|) are the first variations of the 
mean curvature k, which in the literature (see e.g. [3]) are known as the second 
variations of area. They are usually presented in the following equivalent form: 



k 



k= £ - 
2 



ki) 



R(9ab) + k AB k AB + k 2 )r + 2rl 



0(e 2 



(96) 
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